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1 Introduction 



One of the early milestones of space physics was Stormer's theoretical analysis 
of charged particle motion in a purely magnetic dipole field [1,2]. This seminal 
study provided the basic physical framework that led to the understanding of 
the radiation belts surrounding the Earth and other magnetized planets. The 
radiation belts are now known to be composed of individual ions and electrons 
whose motion is often well described by magnetic forces alone. These classical 
results are also relevant to the dynamics of charged dust grains in planetary 
magnetosphercs. However, the much smaller chargc-to-mass ratios produce 
a more complex dynamics, as planetary gravity and the corotational electric 
field must also be taken into account [3-6] . 

In a series of recent papers [7-10] equilibrium and stability conditions 
were derived for charged dust grains orbiting about Saturn. These orbits can 
be highly non-Keplerian and include both positively and negatively charged 
grains, in prograde or retrograde orbits. The first article was restricted to equa- 
torial orbits, while the second treated nonequatorial "halo" orbits, i.e. orbits 
which do not cross the equatorial plane. Both assumed Keplerian gravity, an 
ideal aligned and centered magnetic dipole rotating with the planet, and con- 
comitant corotational electric field. The third paper dealt with the effects of 
planetary oblateness (J2), magnetic quadrupole field, and radiation pressure. 
While the first two forces were found to have a neghgible effect on particle 
confinement, the effects of radiation pressure could be large for distant orbits. 
Interestingly, J2 and radiation pressure can act synergistically to select out 
one-micron grains in the E-Ring [11]. The final paper in this series allowed the 
surface potential of a grain, and hence its charge, to adjust to local photoelec- 
tric and magnetospheric charging currents. It was concluded that stable halo 
orbits were mostly likely to be composed of rather small (~ lOOnm) positively 
charged grains in retrograde orbits. A dust grain "road map" was drawn for 
the Cassini spacecraft now en route to Saturn, showing where to expect dust 
grains of a given composition and radius. 

This paper presents a more comprehensive treatment of dust grain dynam- 
ics, but under the simplified assumptions of Keplerian gravity, pure dipole 
magnetic field, and no radiation pressure. Some of the results were already 
presented in the letters [7,8], here we fill in the details of the necessary calcula- 
tion and also present new results. Our goals are a mathematically rigorous yet 
simplified derivation of equilibrium and stability conditions which highlights 
the relative importance of the several different forces acting on an individual 
grain. 

As is well known, there are no stable equilibrium circular orbits for the pure 
Stormer problem of charged particle motion in a pure dipole magnetic field. It 
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is the addition of planetary gravity and spin that gives rise to stable families 
of equatorial and nonequatorial orbits. We begin with a general discussion of 
charged particle motion in axisymmetric geometry, which is then specialized 
to the motion of charged grains in a planetary magnet osphere. Equilibrium 
conditions are derived first for equatorial orbits, then for halo orbits. Next 
we take up the issue of stability for each family of equilibrium orbits. Results 
are presented for four distinct problems: the Classical Stormer Problem (CSP) 
in which a charged particle moves in a pure dipole magnetic field, the Rota- 
tional Stormer Problem (RSP), with the electric field due to planetary rotation 
included, the Gravitational Stormer Problem (GSP), with Keplerian gravity 
included but not the corotational electric field, and the full system (RGSP) 
including both fields. For each case one must also consider each charge sign in 
prograde or retrograde orbits. Our results may be summarized as follows: 

CSP: As is well known, no stable circular orbits, equatorial or nonequatorial 
exist. However, under adiabatic conditions important families of guiding 

center orbits confined to a potential trough called the Thalweg exist. 
Such trajectories lie outside the scope of the present paper. 

RSP: Stable equatorial equilibria exist for both charge signs. There are no 
halo orbits. 

GSP: Stable equatorial equilibria exist for both charge signs. Positive halos 
are retrograde and negative halos are prograde. Both types are stable 
wherever they exist. 

RGSP: Stable equatorial equilibria exist for both charge signs. There is a 
range of positive charge-to-mass ratios without stable equatorial equi- 
libria. Negative halos are prograde, while positive halos can be pro- or 
retrograde. For stability the frequency must be sufficiently different from 
twice the rotation rate of the planet. 

Therefore halo orbits do exist with and without the corotational electric field. 
However, the corotational electric field is required in order to sustain stable 
positive prograde halos. 

2 Charged particle dynamics in axisymmetric 
geometry 

The equations of motion of a particle of mass m and charge g in M^, r = 

(a;, y, ^)*, are 

mr = -B X r - Vt/(r) , 
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where the potential U (r) generates the forces of gravity and perhaps corota- 
tional electric field. Denote by a rotation around the z axis and assume that 
the magnetic field B and the potential U are symmetric with respect to this 
rotation: 

B{Rr) = RB{r), VU{Rr) = RVU{r) . 

In particular this is true for the field B = V x A of a centered magnetic dipole 
of strength Ai and dipole axis the z axis, for which the vector potential is, in 
the Coulomb gauge, 

A(r) ^M{y,-x, Of/r^, = + + 

The equations of motion can be transformed to a rotating coordinate sys- 
tem using a rotation matrix R corresponding to the angular velocity Q — 
(0, 0, ny which rotates around the 2;-axis with with angular speed fl. For 
given angular velocity Q the z axis is chosen in the same direction, so that Q 
is positive. Note that the magnetic moment A4 can be positive or negative. 
Direct differentiation then gives 

r = i?q, r = i?(q + Qxq), r = i?(q + 2Q x q + Q x (Q x q)) . 

so that 

mq = (-B - 2mO) x q - x x q) + -B x x q) - V[/(q) , 
c c 

where —2mQ x q is the Coriolis force and —mfl x (f2 x q) is the centrifugal 
force. Following the standard argument the term x (f2 x q) is not present 
in a frame rotating with the planet, because there is no additional electric 
field. Therefore there must be the same term in Vf/ in order to cancel it. In 
particular we model the situation inside the corotating plasma, for which this 
assumptions is reasonable. Transforming back to the rest frame this addition 
to the potential gives the corotational electric field observed in the rest frame. 
The corresponding potential is obtained from 

2 2 

E = -B X (Q X r) = 7fiV*, * = ^+y_ ^ 
c r"* 

where 7 = qA4/c. The electric field E is unipolar, i.e. its curl is zero, and 
therefore it is not induced by a changing magnetic field. Moreover, it is per- 
pendicular to the magnetic field, E • B = 0. The divergence of this field is not 
zero; instead we find V ■ E = 2jQ{2z'^ — — y^)/^^i so that there is a space 
charge distribution originating from the rotation of the plasma. 
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In an inertial frame the potential now reads 



U{r) ^-ag— + (7^7Q* , 



where the parameters cr^ and 0"^ serve as markers of the gravitational and 
electric forces in order to track the origin of the various terms after scaling away 
excess parameters. Usually we consider cr^ = (Xp = 1; the classical Stormer 
problem [1] has a^. — (Jg — 0. The case (7^ = and Og — l has also been studied 
by Stormer [1,2]; we shall refer to it as the "Gravitational Stormer Problem." 
While it is simpler than our case, it will turn out that the most important 
physical effects can already be seen in this subcase. Another interesting special 
case is = and — 1. This takes into account the effect of the electric field 
in the rotating plasma, but the particle is still massless. We call this case the 
"Rotational Stormer Problem." It will turn out that in this case halo orbits 
do not exist. 

The Hamiltonian of the above equations of motion is 



Owing to the symmetry of the problem the z axis is an invariant set. Because 
it is singular in coordinates adapted to the symmetry of the problem it is 
best analyzed in cartesian coordinates. The magnetic field is parallel to this 
axis, so that for motion on this axis there is no Lorentz force. For initial 
conditions in the set x — y — p^ — Py — the equations of motion in cartesian 
coordinates show that the derivatives of x,y,px,Py are zero; hence it is an 
invariant set, on which the Hamiltonian is purely gravitational, Hz = pl/2m — 
Hm/z. Depending on the initial conditions in the invariant set a particle either 
collides with the planet or escapes to infinity; there are no stationary points 
on this axis. 

In cylindrical coordinates (p = ^/x^~+^, (f) — arctan(y, x),z) H becomes 



with the dipole stream function ^! = p^/r^. For systems with symme- 
try a stream function ^ independent of (f) can always be introduced. The 
corresponding vector potential is (sin0, — cos0, 0)*^/p, the magnetic field is 
{x^'z.y'^z, -p'^pY/p^, and the electric field is Vl{x^! p,y^! p, p'^ / P- 
We may distinguish three types of constants in the problem: 

• Parameters describing the planet's mass ji — GM and spin rate Q,. They 
are the most fixed parameters. 
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• Parameters describing the dust particle's mass m and charge, measured 
by 7 = qM. / c. 

• The angular momentum and total energy h = H are constants of 
the motion determined by the initial conditions. Fixing both h and 
defines a region of possible motion in configuration space. 

We now introduce a convenient scaling to reduce the number of parame- 
ters. Time is measured by the inverse frequency of the planetary spin rate Q. 
Distances are measured in terms of the radius of the Keplerian synchronous 
orbit 

while mass is measured in units of the particle mass m. The scaled Hamiltonian 
is then 

H = ^{pi+Pi+[-.-si,] 1-^ + ^?^, (3) 

where the variables with hat are measured in the new scale. From now on we 
drop all hats. The essential dimensionless parameters are 

R5 . . il'y q M 
p — P(j) — and = 




7 m/i m c GM ul 

where Uc = qBo/mc is the cyclotron frequency, Bq is the planetary magnetic 
field on the equator, ujk = \/ GM / R^ is the Kepler frequency, with Rg the 
planetary radius, and the parameter p is just the angular momentum mea- 
sured in the new units. The single parameter for the dust grain is 5, which is 
essentially the charge-to-mass ratio. Recall that the z axis is oriented so that 
f2 > 0. In the following we will loosely talk about positive/negative charge 
when we mean positive/negative 6. This correspondence is correct if the mag- 
netic dipole moment Ai is positive, i.e. the spin and the field are aligned. This 
is true for Saturn, the main application that we have in mind. Our results are 
valid in both cases. 



3 Equilibria 

3.1 Equatorial Orbits 

Here we shall find it advantageous to work in spherical coordinates (r = 
— axccos{z / r) , (p) , rather than the cylindrical coordinates of Ref 
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[7,8]. The Hamiltonian becomes 
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H = ^[pi + 'f,)+Ue., (4) 

where the effective potential Ues is the part of the Hamiltonian independent 
of the non-constant momenta Pr and pg: 

{pr- 5 sin^ef ag-ar5sin'^0 

Ues{r, e,p) = — 2 . 5 

2r^ sm r 

The equations of motion are then 

r = pr, 9 =pe, (j) = dpUcs, (6) 
Pr = -drllee, pe = -deUcs, p<j> = . (7) 

In order to facilitate the calculation of the partial derivatives of C/eflf we intro- 
duce the frequency 

^^r,e) = 4> = d,U,s = --■ (8) 

In analyzing circular equilibrium orbits it is preferable to employ uj rather 
than p^ as parameter, as it is the sign of uo that determines whether the orbit 
is rotating in the same direction as the planet (prograde) or opposite to it 
(retrograde). Recall that in the scaled variables frequencies are measured in 
terms of Jl; hence u — 1 means synchronous motion. Prom now on we will 
eliminate p in favour of a; = a;(r, ^) in the potential to get 

TT 12 2-2/) 0-9 , o-r(^sin^6l 

C/eff = -UJ r sm 6 -\ . (9) 

2 r r 

It is important to notice that when calculating derivatives of C/efr with respect 
to r and 6 we have to treat a; as a function of r and 9. The derivatives of uj 
(re-expressed in terms of uj) are 

drUJ = -2 h —r 

r 

dguj — —2 cot 9 iuj + ^ 



so that 



drills = -u'^rs\n^9+^{5{u-ar)sn\^9^ag) (10) 

deU^s = -- cos e sin e(u;V + 2u;(5- 2(7^(5). (11) 
r 
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If all partial derivatives with respect to r, 6, and p are zero there is no 
motion at all. The latter is just u, which we set to zero. Then (p!TD requires 
6 = tt/2. The other possibility = 0, vr is the coordinate singularity, and it has 
already been treated using Hz- The uncharged case 5 = is not of interest here. 
For 9 = ii/2 and u = ([T0|) reduces to ag = arS. Therefore particles at rest can 
occur anywhere in the equatorial plane but only when 6 = (Jg/cJr- Therefore 
5 = 1 can be considered as the case where electrical and gravitational forces 
are balanced. In the classical Stormer problem a,. = o"g = both equations 
are automatically satisfied once = 0. In this case we can place a particle at 
rest anywhere in space; since there is no potential there are no forces if there 
is no motion. In any case it is true that for these trivial solutions at rest the 
angular momentum p is nonzero; from uo = Q and (||) we find p = S/r. The 
energy is zero for all these equilibrium points. 

The system has a discrete symmetry: the equations of motion are invariant 
under the map {0,pg) — > {it — 6, —pe)- The set that is invariant under this 
map is the equatorial plane with no transverse momentum, {6,pg) = (7r/2,0). 
As always this is also an invariant set for the dynamics. The physical reason 
for this is that B(r) is parallel to the z axis if r = {x, y, 0), so that for motion 
within the equatorial plane there is a Lorentz force, however, with direction 
in the plane. Moreover, E(r) is perpendicular to B(r), and therefore E(r) 
lies in the equatorial plane if r = {x,y,0). The Hamiltonian restricted to the 
equatorial plane reads 



Hxy 




ag — a J, 5 



This is an integrable system with one effective degree of freedom that can be 
solved in terms of elliptic functions. The effective potential in the equatorial 
plane is 

2 \r r^J r 2 r 

The minima rg of Uxy{T) correspond to circular orbits in the equatorial plane 
because the right hand sides of all the equations of motion except cj) = uj are 
zero. The calculation of critical points of Uxy{r) leads to the solution of a cubic 
polynomial in r given by r'^drUxy If instead we eliminate p in favour of uj we 
obtain a much simpler polynomial 

P{r,uj) = uj'^r^ - uj5 + ar5 - Ug. (12) 

Solving P = for r yields a generalization of Kepler's 3rd law for equatorial 
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Figure 1: Existence of equatorial orbits 



orbits: 

rJ^y^^-All^ZjlA. (13) 

which for (5 — reduces to the ordinary Kepler's law. The corresponding 
angular momentum p can be calculated from and is Pe = ujr"^ — 5/rf.. The 
radius is positive if 

<ar-'^ and 5 < 0, or (14) 



io>(Xr-^ and 5>0. (15) 



For negative 6 there are always pro- and retrograde orbits. For positive 6 
this is only true for 6 < Ogjor, while for 5 > Ogjor all orbits are prograde. 
In Fig. |l] the possible combinations of lo and 5 for which circular equatorial 
orbits exist are shaded grey. The horizontal asymptotes have io = Or- The 
hyperboloidal boundaries and the limit u; oo correspond to zero radius. In 
the Stormer case = = all equatorial circular orbits of negatively charged 
particles (5 < 0) are retrograde (cu < 0), while positively charged particles have 
prograde orbits. The gravitational and/or electric field perturbations create 
small regions of the opposite behaviour for some 5. In addition some motions 
for large positive 8 and small positive uj are made impossible by switching on 
the additional fields. 

The equation for equatorial orbits, P = 0, can be solved for 5 in order to 
give b function of u) for given r: 

= . 

uj — 

These curves are shown in Fig. |^ and ^. Note the two straight lines with 
r = 1, which in our scaling is the radius of the synchronous orbit in the Kepler 
problem. The horizontal one corresponds to the synchronous orbits {uo = 1) 
which exist for any 6. Hence the synchronous Kepler orbit is not affected 
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Figure 2: Curves of constant radius r = O.li, i = 1, . . . , 20 for equatorial orbits. 
The thick hues have r = 1. 

by the addition of both fields. This is not true for the three Stormer cases, 
see below. Another prominent feature is the point {S,uj) = (1,0), which is 
intersected by hyperbolas with all r. It corresponds to the equilibrium points 
discussed above. 

Let us briefly discuss the corresponding diagrams for the three Stormer 
cases shown in Fig. ^. In the classical case u and 6 are proportional. Small 
slope means large r. The relation breaks down for the line a; = for which 
all radii are possible. There are only two other types of equatorial orbits: 
positively charged prograde and negatively charged retrograde. In the gravi- 
tational Stormer case the relation between u and S is quite similar, but now 
there are also small regions of the two types of motions: positively charged 
retrograde and negatively charged prograde. In the rotating Stormer case the 
positive retrograde orbits have disappeared again. The main new feature is the 
appearance of orbits of small negative charge with small u and small radius. 



3.2 Halo Orbits 

Our goal is the calculation of periodic orbits that encircle the planet in a 
plane parallel to the equatorial plane but entirely above/below it. Circular 
orbits correspond to critical points of Ugs, i.e. points (ro,^o) which both 
derivatives of f/efr vanish. This is so because at {r,9,Pr,Pe) = ('"o;^0;0,0) the 
right hand sides of Hamilton's equations are zero, except for (p. They are given 
by the minima of t/efr. Their stability will be calculated in the next section. 
Circular orbits are given by the solution of drU = and dell = for 



arbitrary u, see ([TO| , p!T|) . The second equation has the solution 6 = n/2, 
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Figure 3: Curves of constant radius for equatorial orbits of the a) classical = 
ag = 0), b) gravitational (ur = 0, cr^ = 1), and c) rotational {a^ = l,(Tg = 0) Stormer 
problem. Thick lines have r = 1. 



which gives the equatorial orbits we already analyzed. Also the solutions with 
6 = 0,n have already been described using the reduced Hamiltonian Hz- The 
remaining solutions are given hj Q = with 

Q{r,uj) =uj^r^ + 2uj5 -2ar5 , (16) 

which describe the nonequatorial (or halo) orbits. The equation Q = can be 
solved for r^, which can then be eliminated from ([TT|) resulting in an angular 
equation A = with 

A{e,uj) = ag + 35{uj-ar)sm^e. (17) 

The functions Q and A completely describe the halo orbits. In particular 
these equations can easily be solved for r and 9, so that all circular orbits are 
obtained in parametric form, with u; as a parameter. Explicitly we find 

thH' = 26^^ (18) 
sm'e.iu) = . (19) 

grains [4]. The second equation clearly shows that without the gravitational 
forces {(7g = 0) there are no halo orbits. In particular in the classical Stormer 
problem there are no halo orbits (except for the trivial equilibrium points 
discussed above). Adding the electric field alone there still are no halo orbits. 
In both equations it is obvious that the electric field merely shifts the frequency 
uj of circular halo orbits. We conclude that the electric field is not essential for 
the existence of halo orbits. 
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Figure 4: Existence of lialo orbits 



For halo orbits the essential condition for their existence is that sin^ 6' < 1, 
which implies 

> f^r - S and 6 <0, or (20) 

^ < f^r - S and 6>0. (21) 

These conditions automatically imply that the corresponding radius is posi- 
tive. Note that the ordering of the u inequalities is reversed compared with the 
equatorial case (14,15). Hence for negative charge only prograde orbits exist 
while for small positive charge only retrograde orbits exist. Only if cTj, ^ and 
for 6 > ag/ScTr can both types of orbits exist. The electric field does make a 
difference for the existence of synchronous halo orbits: If cTj. = 1 then u = 1 is 
impossible for finite 6. If cr^ = synchronous halo orbits exist for 6 < —ag/3. 
Note that u = ar (in particular a synchronous orbit if cxr = 1) is impossible for 
finite S. Without the electric field there do exist synchronous halo orbits with 
negative charge. In order for (g/c)B x r to balance gravitation for positive 
charge we have to reverse r, hence u = —1. 

To get an overview of all possible halo orbits we plot curves of constant r 
and 6 in the 6-uj plane. The family of curves of constant radius is given by 



6 



2{ar — uj) 

The family of curves of constant azimuth 6 is given by 



S- 



3((7.r — uj) sin 6 

Both families are shown in Fig. ^ for a.,. = cXg = 1 and in Fig. ^ for the Stormer 
problem with gravitation. 

The regions of existence in the S-u-plane for equatorial and halo orbits 
only overlap in a small region bounded by hyperbolae. They are in a sense 
connected at the hyperboloidal boundary of the halo orbits, because in the 
next section we will find that this line marks a pitchfork bifurcation of an 
equatorial orbit changing its stability and creating halo orbits. 
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Figure 5: Curves of constant r and Q for halo orbits. Thick lines have r = 1. 
r = 0.2i, sin (9 = O.li, z = 1, . . . , 10 
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Figure 6: Curves of constant r and for halo orbits in the gravitational Stormer 
problem. Thick lines have r = 1. r = 0.2i, sin0 = O.li, « = 1, . . . , 10 
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4 Stability 



In Ref [7-8] explicit stability boundaries for both equatorial and halo orbits 
were calculated. Here we obtain these boundaries more directly using the 
fact that all circular orbits may be parameterized hj u. A circular orbit is 
stable if it corresponds to a local minimum of U, for which we need the second 
derivatives of U. Using the chain rule gives 

fftU = (3^^-25^^ + C^sin^- 



r 



4.1 Equatorial Orbits 

For equatorial orbits we insert 9 = ti /2 and r = re{uj) as given above. The 
mixed derivative vanishes and the other two are 



2 \Og + 5{uj - (T,.; ' 

dlU\e = ag + 35{uj — ar) . 



The radial derivative diverges for those uj that correspond to r = 0. It vanishes 
for UJ = 0; however, the corresponding radius is not finite (except for 6 = 
(7g/(Jri the case of equilibrium points). There are two nontrivial factors that 
correspond to tangent bifurcations of equatorial orbits with 

± 1 ± 73 / ag 

The vanishing of the second 6 derivative indicates the loss of transverse sta- 
bility. Because of reflection symmetry 9 —9 this results in a pitchfork 
bifurcation with 

UpF — Cr-r — 777- 

3o 

All three curves are hyperbolas in the space of {6, u) shown in Fig. |^. The 
intersection of the curves upp and occurs at 

± ^ gg 5 ± 2^3 
ar. 3 
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Using the above formulas for the second derivatives it is easy to check that 
stabihty only holds in the following a; ranges (compare Fig. |^). 

6 < : < u < ojpF 
< 5 < 5^ : ujpp < oj < uj~ 
5 > 5'^ : uJpF < OJ < oj'^ . 

For the first two ranges = is (partially) included, which means that the 
corresponding family of orbits exists for arbitrary large radius. The radius r^. 
as given by ([TsD as a function of is a monotone function for most of these 
orbits except for unstable orbits with 5 > 1. This gives the additional curve 
in the diagram. The fact that the radius is not monotonic can already be seen 
in Fig. 1^, where the turning points are marked by a gray line. 

For comparison we also show the diagrams in the two nonclassical Stormer 
cases, see Fig. ^. In the classical Stormer problem with = = the deriva- 
tive 9^f/eff|e is always negative, hence there are no stable equatorial orbits. In 
the purely gravitational case stable orbits exist between the pitchfork curve 
and the tangent bifurcation corresponding to uj~ . The curve cu^ always is out- 
side the region of existence. In the rotational Stormer problem the pitchfork 
curves and the existence curves coincide at = 1. The stable region for 5 < 
is between uj~ and 1, while for 5 > it is between 1 and . Comparing these 
pictures with Fig. |^ one clearly observes that for large \5\ (i.e. small mass) 
the systems behaves like the rotating case, while for small \5\ (i.e. large mass) 
the behaviour is dominated by gravity and looks like the gravitational case. 
In the gravitational case there exist stable pro- and retrograde orbits for any 
5; the system is symmetric with respect to change of sign of uj and 5. In the 
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Figure 8: Regions of stability for equatorial orbits of the a) gravitational, and b) 
rotating Stormer problem are dark. Overlay Fig. ^ 

rotational case this symmetry is broken and both pro- and retrograde orbits 
only exist for negative charge. Trying to interpolate between the large and 
small 6 behaviour the full system creates an interval [6~ , 5^] of charge-to- mass 
ratios for which no stable equatorial circular orbits exist. 



4.2 Nonequatorial Orbits 

In Ref [8] the stability of halo orbits was analyzed by examining the zeros of 
a quintic in p. Here we obtain similar results more directly. For nonequatorial 
orbits the calculation is essentially the same as for their equatorial cousins, 
except that all three second derivatives are nonzero and we have to calculate 
the determinant and the trace of the Hessian, which turn out to be 



detD^Ueslh 



2u?{u? - 4uj(Tr + cr^)(ag + 3S{uj - cr,)) 
35(cj- 1)3 



-0", 



Or 



The determinant vanishes in three cases: 1) For = 0, which again for 5 7^ 1 
does not correspond to finite orbits. Since lo appears squared there is no change 
in stability when orbits go through infinity. 2) In the case that the last factor 
is zero, which reproduces the condition a; = ujpp- 3) There are two new critical 
10 given by the remaining factor as 



These frequency values correspond to tangent bifurcations of halo orbits. The 
corresponding horizontal lines are also shown in Fig. |^. A pair of stable and 
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unstable orbits is created for this frequency. The hues only extend up to the 
intersection with the pitch fork curve, which occurs at 



6 



~ ar 6 



For 6 < the upper curve extends up to 6^ , the lower curve is valid for 6 > 
and extends from 6'^ to infinity. For 6 < halo orbits have to be above the 
pitchfork line. Inserting into the invariants of the Hessian we find that they 
are stable if they are above u = and unstable otherwise. In the unstable 
family there occurs a maximum in radius at = 2. Otherwise the radius 
is a monotonous function of u. For 6 > stability is reversed: orbits exist 
below the pitchfork line and are stable if u < . At passage through u = 
the radius goes to infinity, so that (sufficiently) positively charged retrograde 
halo orbits exist for all large radii. Hence we obtain the following u ranges of 
existence of stable halo orbits: 



(5 < 5_ 
6-<6 <0 

6>6^ 



UJ > UJpF 
UJ < UOpF 



A simple way to characterize Fig. ^ is to say that halo orbits with frequencies 
too close to synchronous are unstable. However, the range of unstable frequen- 
cies is centered around orbits with uj = 2, i.e. orbits that go around twice for 
one revolution of the planet. Halo orbits with frequencies further away than 
-\/3 from this are stable. Note that the equatorial orbits behave approximately 
in the opposite way. For them only orbits with small uj are stable (except for 
small 6). This can also be interpreted in terms of the pitchfork bifurcation: 
once equatorial orbits become too fast, they become unstable and create stable 
nonequatorial orbits. The corresponding picture for the gravitational Stormer 
problem is not shown, because it is trivial: in this case every halo orbit is 
stable, as can be seen from the above expressions for determinant and trace of 
the Hessian. 



5 Stable Halo Orbits in Space 

Considering Fig. ^ we see that the curves of constant r and 9 transversely 
intersect each other in the regions of existence. This means that the transfor- 
mation from (r, 9) to (5, uo) is invertible, which we will now show. Instead of 
transforming to spherical coordinates (r, 9) we directly transform to cylindrical 
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Figure 9: Regions of stable halo orbits are dark. Overlay Fig. |5| 

coordinates. Equations (|T8|) and (|T9|) can be considered as a transformation 
from {z,p) = {r cos 9, r sin 9) to For each of the four types of orbits 

distinguished by pro/retrograde and positive/negative charge this is a global 
transformation because the Jacobian is 



det 



diz,p) _ 2 
d{6,uj) 96uj^r sin 9 cos 9 



which is only singular when lj, 6, r, or sin 29 is zero. We already know that the 
latter two are only zero at the boundaries of the valid region in {S, ut) space. 
The inverse of the transformation is given by 



1 



3 sin^ 9 



1 

3 (ar 



1 



— io) sin 9 



(22) 
(23) 



The first is a generalization of Kepler's 3rd law for halo orbits, which surpris- 
ingly is independent of the electric field. Compared to the usual law it has 
effective frequency ^^/3/2uj sin 9. From the second equation the corresponding 
charge to mass ratio can be calculated. These two equations give a precise pre- 
diction about what dust particles with what velocities should be observable at 
a given nonequatorial position (r, 9). 

To get an idea about what particles to expect at what position we now plot 
the curves of constant uj and 6 on the p-z plane. The simplest way to do this is 
to use (|TB|,|T^ to generate a parametric form of these curves in the p-z plane: 

{p,z) = {rh sin 9h,rh cos 9h). 
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Figure 10: (5, grid in {p,z) space for halo orbits in the gravitational Stormer 
problem, (top left to right ±uj = -2 - ^/2>, -1, -0.5, -2 + ^/S, -0.15, -0.1, -0.08, 
bottom ±5 = 0.1, 0.5, 1, 2, 3, 4, 5, 6) 
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Figure 12: ((5, w) grid in (p,z) space for a) prograde positive (w = 2 — 
\/3, 0.15, 0.1, 0.08, left top to bottom S = 1.5,1,0.8,0.6,0.5,(1 + \/3)/6, 0.4, 0.37), 
and b) prograde negative {lo = 15, 8, 5, 2 + -y/S, right top to bottom 6 = 
-50, -25, -10, -3, -1, -0.3, (1 - \/3)/6, -0.06, -0.03) halo orbits. 



In the gravitational Stormer case the formulas are 
25 ^ 



ip,z)=(--\ (v/-1/(3M,V1 + 1/(3M)- 



6 and u have to be restricted to the range of existence, respectively stability, 
which is the same in the present case. The resulting diagram is shown in 



Fig. [I^. Because there is no coupling to the rotation, prograde and retrograde 
orbits are the same up to the sign of u. 

For the full system there is an additional region of prograde orbits with 



5 > 0, see Fig. |T2|. In this case only {5,u) values from the stable regions are 
taken to draw the grids. This is the reason for the cutoffs in the prograde case. 
Note that all the retrograde halo orbits are stable, and therefore for every point 
in the (p, z) plane there exists a unique halo orbit. 

Note that all four figures share the same set of lines uj = const. This is a 
result of the fact that the generalized 3rd Kepler's law (^) is independent of 
ar, S, and independent of the sign of u. Using the transformation from (r, 6) 
to (p, z) we can convert this equation into 

2 9 f 2 



which is the explicit form of the curves uj = const in all four figures. A similar 
explicit form for the curves 6 = const can only be obtained with cXr = 0. Then 
we find z'^ = (— v^5p)^/^ — p^. With cr^ = 1 the curves in (p, z) are described 
by a polynomial of degree 5 in z"^ and p^, so that the above parametric form 
is the most convenient representation. 
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Our most important conclusion is that the dependence on uj is the same 
with or without the electric field. The distribution of grain sizes as given by 
the curves of constant b is significantly changed. The changes are fairly small 
for stable retrograde positive orbits. In both cases they exist at any point in 
space. The prograde negative orbits need quite high angular velocity and only 
survive close to the z axis. Stable prograde positive halo orbits do not exist 
at all without an electric field. With the field they need to have a minimal 
distance of a little more then twice the synchronous radius, and b must be 
around (1 + v^)/6 in order to be able to be close to the planet. It follows that 
retrograde positive orbits are the most likely candidates for halo orbits. 

6 Discussion 

We have calculated explicit equilibrium and stability conditions for arbitrary 
circular orbits in an axisymmetric combination of gravitational, magnetic and 
corotational electric fields. The equilibrium and stability boundaries were con- 
veniently parametrized by the charge-to-mass ratio b and the orbital frequency 
uj. The individual effects of planetary gravitational field, magnetic field and 
corotational electric field on the existence and stability of charged dust grain 
orbits were elucidated. 

Our principal result is that halo orbits cannot exist without inclusion of 
gravitational forces. Without the corotational electric field all halo orbits are 
stable. The distribution of orbital frequencies of stable halo orbits in space is 
the same with- and without the corotational electric field, which is the content 
of a generalized Kepler's 3rd law (P^). The inclusion of the corotational electric 
field alone does not give halo orbits at all. Adding it to the gravitational 
field does not have a strong effect on positive retrograde orbits, which are 
still all stable. It destabilizes negative prograde orbits with small frequencies. 
Adding the corotational electric field has a surprisingly strong effect on the 
character of both equatorial and nonequatorial (halo) orbits. In particular, 
prograde positively charged halos require a corotational electric field for their 
very existence. 

For halo orbits lying several Saturn radii above the equatorial plane the 
typical surface potential of a dust grain is expected to be around +5\^, due to 
the low plasma density there and resultant dominant photoelectric charging. 
If stable retrograde grains are present, even the very small grains predicted by 
our theory should be detected by the CDA experiment on board the Cassini 
orbiter due to arrive in 2004. 
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